Low- Temperature Thermal Conductivity of Superconductors With Gap Nodes 
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We report a detailed analytic and numerical study of electronic thermal conductivity in d-wave 
superconductors. We compare theory of the cross over at low temperatures from T-dependence to 
T 3 -dependence for increasing temperature with recent experiments on YBa2Cu307 in zero magnetic 
field for T 6 [0.04^,0.4^] by Hill et al, Phys. Rev. Lett. 92, 027001 (2004). Transport theory, 
including impurity scattering and inelastic scattering within strong coupling superconductivity, can 
consistently fit the temperature dependence of the data in the lower half of the temperature regime. 
We discuss the conditions under which we expect power-law dependences over wide temperature 
intervals. 



Low-temperature measurements of electronic trans- 
port properties have given a wealth of information about 
nodal quasiparticles deep in the superconducting phase 
of unconventional superconductors [ll|2|, |^|j,|5|,|al2lLi 
IE El EH 113 El ■ The ultra-low temperature regime is 
of great interest because response functions such as the 
zero-frequency charge conductivity, thermal conductiv- 
ity, and thermoelectric response function, are believed to 
be limited by elastic impurity scattering and obey power 
laws | 141 Il5l llfjl Il7j . Thermal conductivity is of particu- 
lar importance in this context, since theory predicts uni- 
versality in the sense that the low-temperature asymp- 
totic does not depend on the properties of the impurity 
potential pjj- F° r higher temperatures, of the order of 
the low-energy impurity band width 7 and higher, ther- 
mal conductivity is not universal and is sensitive to the 
details of impurity scattering [l6j |. 

The universal character of the T — > thermal con- 
ductance was studied experimentally in great detail for 
several different cuprate materials for a wide range of 
doping levels 0, 0, H H E3- Measurements of 
how the universal limit is approached as temperature is 
lowered was only recently reported by Hill et al [l2( for 
YBa2Cu307. After subtracting a phonon contribution 
oc T 3 , they found an electronic contribution of the form 
Kel /T = ( Ko /T)(l+fcT 2 ), with k /T = 0.16mW/(K 2 cm) 
and the coefficient k = 19.2K~ 2 . This form is consistent 
with a Sommerfeld-type expansion of « e ;(T) f |l6j and 
below). However, the large value of k implies an unre- 
alistically clean sample if impurity scattering is in the 
unitary limit. On the other hand, for Born limit scat- 
tering, the concentration required to fit the large slope 
is unrealistically large and should lead to a severe reduc- 
tion of the critical temperature which is not seen. Their 
plausible conclusion was therefore that scattering might 
be in between the unitary and Born limits. 

In this paper we present results for thermal conduc- 
tivity including effects of elastic impurity scattering and 
inelastic electron-boson scattering. We focus on how 
the universal limit is approached and make a detailed 



comparison with the experimental results of Hill et al 
[l2| . Our results can be summarized as follows. (1) The 
main effect of inelastic scattering is an effective mass- 
dependence of the universal (with respect to impurity 
scattering) T — > thermal conductance |l9|. The coef- 
ficient k is however effective mass independent. (2) We 
can fit the data |l2j in the lower half of the T— interval 
with a best fit for a phase shift slightly below 90°. (3) 
The discrepancy at higher T is related to the split of 
the low-energy impurity band, which makes it hard to 
observe power laws over wide T— intervals in ultra-clean 
materials. 

The thermal conductivity k(T) of a superconductor 
has an electronic contribution K e i(T), a phononic part 
K p h(T), and possibly contributions from other existing 
excitations. We will only consider n e i(T). To quan- 
tify the effects of inelastic scattering on the same footing 
as elastic impurity scattering we make a self-consistent 
strong-coupling calculation using quasi-classical Eliash- 
berg theory [20| including impurity scattering in the 
t— matrix approximation. We consider a two-spectra 
model for electron-boson interactions that generate both 
d-wave superconductivity and an inelastic scattering self 
energy. One spectrum gives coupling in the s-wave and 
d-wave channels, but is attractive only in the d-wave. 
Hence, this part is responsible for pairing but also con- 
tributes to inelastic scattering. The second spectrum 
couples only in the s-wave channel (repulsively) and gives 
an incoherent background scattering. Even though the 
model is general we choose parameters that generates a 
d-wave superconductor with properties close to that of a 
weakly coupled superconductor, with a spectroscopic gap 
A = 2.59 T c (T c = 0.13w mode ) j^. The temperature de- 
pendence K e i(T) is displayed in Fig.^for several impurity 
concentrations rij mp (ro = n\ mp /TrAff) and near-resonant 
impurity scattering (phase shift 5 = 89.5°). We show 
in panels (a) and (b) the electron-boson coupling func- 
tions, x 2 /(u;), and the resulting density of states, Ns(e). 
The imaginary part of the inelastic self energy gives a 
T— dependent contribution to the scattering rate , Y- m {T) 
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FIG. 1: Temperature dependence of the electronic contribution to 
the thermal conductivity for elastic scattering rates as indicated, 
for a phase shift <5o = 89.5°. Insert (a): the two electron-boson 
coupling spectra used in the calculations. Insert (b): single-particle 
density of states in the superconducting state at T = 0. 



22, 23, 24], while the contribution from elastic impurity- 
scattering, r = Tq sin 2 Jo, is T— independent. The rela- 
tive importance of the two contributions leads to a cross- 
over temperature T* . Impurity scattering dominates in 
the low— T region, while inelastic scattering dominates in 
the high— T region. The cross-over is reflected in K e i(T) 
as the peak at T ~ T* in Fig. □ H3 

Although the imaginary part of the inelastic self- 
energy is small for T -C T*, the real part remains 
important. By using general symmetry relations we can 
write down the renormalization of the energy from both 
elastic and inelastic scattering as 



1 

-Tr 
2 



T 3 t R (e)} = i 1 +(l + \)e + 0[e 2 ] . (1) 



Here is 7 solely due to impurity scattering. The slope of 
the real part of e R (e) has contributions from both impu- 
rity scattering and inelastic scattering, A = Ai mp + Aj n . 
Because 5RS^(e) is linear in energy up to the mode- 
energy, w mo de > Ao 3> T*, an effective theory can be 
made for T <C T* by dividing the transport equation by 
1 + Ai n evaluated in the superconducting phase. The in- 
elastic part is then interpreted as an effective mass via 
m*/m = 1 + Aj n - The contribution from impurity scat- 
tering on the other hand, has a non-trivial energy depen- 
dence for e < Ao, and can not be divided out. Within 
the effective theory, we use Eq. Q and the method of 
Graf et al to write down the low— T form 



Kel (T) _ n 2 2M f v 2 
T 3 7t/jAo 
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Here are A// and Vf the effective density of states at 
the Fermi level and the effective Fermi velocity, while 



Ao is the spectroscopic gap and /i is the opening rate 
of the gap function at the node. We now find that 
the remaining effect at T <C T* of inelastic scattering 
within strong-coupling superconductivity is a modifica- 
tion of the T — ► asymptotic of the thermal conductance. 
When we write explicitly, A} — > N* f = M9(1 + X;„), and 
Vf — > v*j = Vj/{\ + Ai n ), one factor 1 + X m remains in 
the denominator. Within the bare theory, this result can 
be traced back to that the spectroscopic gap Ao in the 
weak-coupling limit is replaced by the strong coupling 
off-diagonal function W(0) [19j. In our model W(e) is 
a weak function of energy and to a good approximation 
W(0) A (l + Ai n ). This effective mass dependence 
is in sharp contrast to the normal state, in which case 
K %(T) = T(ir 2 /3)J\f f v 2 f /2T is independent of to*, be- 
cause of canceling factors of (1 + A, n ) in Nf,v 2 and T in 
accordance with gen eral strong-coupling results for trans- 
port coefficients |25j . In this sense, transport properties 
of quasiparticles in the low-energy nodal impurity band 
is quite different from quasiparticles in the normal state. 

Effective-mass factors from a and 7 cancel in the coef- 
ficient k = -^-(a/7) 2 in Eq. J5J and it is only dependent 
on properties of the impurity potential. In addition to 
the electron-hole symmetric part above, also the electron- 
hole asymmetric part of the impurity self-energy enters 
|2(j| . Only the imaginary part 
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contributes to the conductivity through a 2 = (1+Ai mp ) 2 + 
2c 2 ^3 ■ The presence and importance of c represents the 
fact that electron-like and hole-like quasiparticles have 
different transport lifetimes above and below the Fermi 
level, except in the strict unitary and Born limits for 
which it vanishes. 

We now discuss the experiment [T3 | in terms of the 
presented effective theory and its low— T expansion in 
Eq. PJ with the parameters {7, Ai mp , c} computed self- 
consistently |l7| . Starting for simplicity with s-wave im- 
purity scattering, we present in Fig. [2Jd) several com- 
binations of impurity concentrations and phase shifts 
that give good fits to the data in the lower half of the 
T— interval. In fact, there is a broad range in parameter 
space for which the slope k = 19.2K~ 2 can be repro- 
duced [blue dash-dotted lines in Fig.|2fa)-(c)]. Similarly, 
we can find many combinations yielding other values of 
k (black and red lines). However, not all combinations 
of parameters produce a consistent fit, in the sense that 
the expansion Eq. (J2J represents over a wide tempera- 
ture range the full numerically computed thermal con- 
ductivity shown in Fig. G{d). The deviation from true 
T 2 -behavior is also seen in the figure: the three curves 
start at low T with the same coefficient, k — 19.2K~ 2 , 
but deviate at higher T. This is also the reason for the 
discrepancy between theory and experiments in the up- 
per half of the T— interval. The part in parameter space 
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FIG. 2: (a)-(c) The low-temperature parameters in Eqs. and 10 for three values of the T 2 -coemcient in Eq. J2J as function of 
scattering phase shift, (d) Three sets of parameters that represents our best fits to the data of Ref. Il2j . The blue shaded region indicates 
the temperature range where a good fit can be made. The dashed curves are from a strong-coupling calculation while the underlying solid 
lines are generated by the effective theory. Curve (II) is the low— T part of the black solid line in Fig.0 In all cases k = 19.2K~ 2 . Note 
that the material appears cleaner in the effective theory as Tq = ro/(l + Ai n ) with To = n\ mp /nAff ■ We normalized K e ;(T)/T to the T — » 
extrapolated value 0.16mW/(K 2 cm). The energy scale is T c = 90K. 



for which the Sommerfeld expansion is applicable can be 
broadly defined as phase shifts for which A; mp is negative. 

Let us discuss in more detail why Eq J2J works 
well only in a restricted low-T interval. Recall that 
the thermal conductivity can be written as n e i(T) = 
J^. deK(e) (—737) , where / is the Fermi function. The 

derivative of the Fermi function is sharply peaked at the 
Fermi level in an interval of order T, which leads us to 
expand all quantities around e = 0, as in Eqs. an d 
©, and do a Sommerfeld expansion K e i(T) — K(0) + 



00 d Zn K{t) 



71=1 



a n T , where the for us relevant coeffi- 



cients a\ — 7r 2 /6 and 02 = 77r 4 /360. However this scheme 
works well only if the function K(e) is slowly varying on 
the scale T. In our case we can write K (e) = e 2 «(e) and 
study the energy dependence of £(e), see Fig.^d). The 
fine structure at low energies for phase shifts below 90°, 
can intuitively be understood by studying the phase shift 
dependence of the self-energy. In a <i-wave superconduc- 
tor with a single impurity scattering strongly (near the 
unitary limit) there is a resonance in the i— matrix at [27| 
e res w ±7rcot<5n/[21n(8/7rcot(5o)] where the sign corre- 
sponds to the electron and hole sectors. The width of 
the resonance is smaller than the resonance energy by a 
factor l/ln(8/7rcot(5o) <C 1. Within transport theory we 
employ an impurity averaging technique [2^] to obtain 
the self-energy for a material with a small concentration 
of randomly distributed impurities. As a result of the 
configuration average, the single impurity resonance is 
broadened into an impurity band of width jb, centered 
around an energy slightly shifted from the single impu- 
rity resonance e res — > £b- The impurity bands are shown 



in Fig. E3a)-(c). For clean systems with a phase shift 
deviating substantially from 90°, the resulting bands are 
split into two parts symmetrically positioned above and 
below the Fermi level. The critical point is when the 
slope at zero energy of the real part of the self energy 
changes sign, i.e. when Aj mp — — •RS^O) turns posi- 
tive. The structures in the self-energy lead to structures 
also in the response functions at low temperatures on the 
scale T ~ e& <C 76, see Fig.^d), and the applicability of 
Eq. (|2J) is drastically reduced in temperature by the new 
constraint T ef,. 

In a more realistic model, impurity scattering is 
anisotropic. The main technical complication in this case 
is the necessity to take care of vertex corrections. They 
are small in the limit T <C 7 |l8j |. but are important for 
higher temperatures. In particular they contribute to the 
coefficient a in Eq. J2J . For definiteness we have included 
p-wave and d-wave scattering channels and considered ar- 
bitrary phase shifts |29j . In general we obtain a pair of 
resonances, one electron-like and one hole-like (symmet- 
rically positioned on each side of the Fermi level), for 
each phase shift. This fact was also noted recently in 
Ref. 301. The main conclusion relevant in the present 
discussion is that the impurity self-energy aquires more 
structures (related to the new resonances) and the range 
of applicability of Eq. J2J is not improved. 

In summary, let us list our results and make a few re- 
marks. (1) Inelastic scattering contributes to the T — » 
asymptotic (i.e. ko oc T), but plays no role for the 
T 3 coefficient (i.e. k). (2) The low-T half of the ex- 
perimental data can be consistently fitted with param- 
eters representing a wide range of impurity concentra- 
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FIG. 3: (a)-(b) The impurity scattering self-energy and (c) the 
density of states for To = 0.01T C and phase shifts ranging from 90° 
(red dashed line) to 80° (thick black line) in steps of two degrees, 
(d) The kernel k{t) of the thermal conductivity contains structures 
at low energies for phase shifts below 90° because of the split im- 
purity bands seen in (a)-(c). Also shown is the derivative of the 
Fermi function which has a peak at the Fermi level of width ~ T. 
Here is T = 0.005T C = 0A5K for T c = 90K. 



tions and phase shifts. The best fit to the large coeffi- 
cient k — \%.2K~ 2 is obtained for a rather clean material 
with a scattering phase shift slightly below the unitary 
limit (So < 90°). (3) However, we find deviations in 
the high— T region. The reason is that the expansion 
Eq. (0) has a limited range of applicability, smaller than 
the temperature range in which the T 3 term was found 
experimentally to dominate. Indeed, in an earlier discus- 
sion of the data ^3] m Ref. [H| , model parameters were 
chosen so that the expansion Eq. was inapplicable. 
In their work the unit-term Eq. © was neglected, which 
changed their results by about a factor two. We con- 
clude that a clear cross over between T— dependence to 
T 3 -dependence, with the respective power laws holding in 
wide T— intervals, is most easily observable in dirty to in- 
termediately clean, but not in ultra-clean, materials. Fi- 
nally, an alternative way of extracting information about 
impurity scattering from transport properties, would be 
to study the thermoelectric effect jl/jj. In this case, the 
non-universality of the T — ► asymptotic can be used. 
This quantity, the analogue of kq, should be easier to 
measure and interpret than the T 3 — part of n e i. 
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